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Abstract 

We obtain possibly valuable information about the orientation-tuning of phase diagram of su- 
perdense nuclear matter at high fermion as well as boson number density but low temperature, 
which is not accessible to relativistic heavy ion collision experiments. Our results resemble those 
proposed before by Alford. Possible observational signatures associated with the theoretically 
proposed states of matter inside compact stars are discussed as well. 
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I. INTRODUCTION 



The only place in the universe where we expect sufficiently high densities and low tem- 
peratures is compact stars, also known as 'neutron stars', since it is often assumed that they 
are made primarily of neutrons (for a recent review, see [1]). A compact star is produced in 
a supernova. As the outer layers of the star are blown off into space, the core collapses into 
a very dense object. In a broader perspective, neutron stars and heavy-ion collisions pro- 
vide access to the phase diagram of matter at extreme densities and temperatures, which is 
basic for understanding the very early Universe and several other astrophysical phenomena. 
These range from nuclear processes on the stellar surface to processes in electron degen- 
erate matter at subnuclear densities to boson condensates and the existence of new states 
of baryonic matter such as color superconducting quark matter at supernuclear densities. 
More than that, according to the strange matter hypothesis strange quark matter could be 
more stable than nuclear matter, in which case neutron stars should be largely composed 
of pure quark matter possibly enveloped in thin nuclear crusts. Neutron stars and white 
dwarfs are in hydrostatic equilibrium, so at each point inside the star gravity is balanced by 
the degenerate particle pressure, as described mathematically by the Tolman-Oppenheimer- 
Volkoff equation [2]. 

It is often stressed that there has never been a more exciting time in the overlapping areas of 
nuclear physics, particle physics, and relativistic astrophysics than today [3]. Neutron stars 
are dense, neutron-packed remnants of massive stars that blew apart in supernova explosions 
[1]. They are typically about twenty kilometers across and spin rapidly, often making several 
hundred rotations per second. Many neutron stars form radio pulsars, emitting radio waves 
that appear from the Earth to pulse on and off like a lighthouse beacon as the star rotates 
at very high speeds. Depending on star mass and rotational frequency, gravity compresses 
the matter in the core regions of pulsars up to more than ten times the density of ordinary 
atomic nuclei, thus providing a high pressure environment in which numerous subatomic 
particle processes compete with each other. 

The most spectacular ones stretch from the generation of hyperons and baryon resonances to 
quark deconfinement to the formation of boson condensates. There are theoretical sugges- 
tions of even more exotic processes inside neutron stars, such as the formation of absolutely 
stable strange quark matter, a configuration of matter more stable than the most stable 
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atomic nucleus, 62 Ni. Instead these objects should be named nucleon stars, since relatively 
isospin symmetric nuclear matter-in equilibrium with condensed K~ mesons-may prevail 
in their interiors [3], hyperon stars if hyperons (E,A,S, possibly in equilibrium with the A 
resonance) become populated in addition to the nucleons, quark hybrid stars if the highly 
compressed matter in the centers of neutron stars were to transform into u, d, s quark mat- 
ter, or strange stars if strange quark matter were to be more stable than nuclear matter. Of 
course, at present one does not know from experiment at what density the expected phase 
transition to quark matter occurs. Neither do lattice Quantum ChromoDynamical (QCD) 
simulations provide a conclusive guide yet. From simple geometrical considerations it follows 
that, for a characteristic nucleon radius of ~ lfm, nuclei begin to touch each other at 
densities of ~ (47rr^/3)) _1 ~ 0.24 fm 3 , which is less than twice the baryon number density of 
ordinary nuclear matter, p = 0.16 fm 3 (energy density e = 140 MeV/fm 3 ). Depending on 
the rotational frequency and stellar mass, such densities are easily surpassed in the cores of 
neutron stars so gravity may have broken up the neutrons (n) and protons (p) in the centers 
of neutron stars into their constituents. The phase diagram of quark matter, expected to be 
in a color superconducting phase, is very complex [3-5]. At asymptotic densities the ground 
state of QCD with a vanishing strange quark mass is the color-flavor locked (CFL) phase. 
This phase is electrically charge neutral without any need for electrons for a significant range 
of chemical potentials and strange quark masses. 

Quite recently McLaughlin et al. [6] searched for radio sources that vary on much shorter 
timescales. They found eleven objects characterized by single, dispersed bursts having du- 
rations between 2 and 30 ms. The average time intervals between bursts range from 4 min 
to 3 h with radio emission typically detectable for < 1 s per day. From an analysis of the 
burst arrival times, they have identified periodicities in the range 0.4-7 s for ten of the eleven 
sources, suggesting origins in rotating neutron stars. Meanwhile as all pulsars from which 
giant pulses have been detected appear to have high values of magnetic field strength at 
their light cylinder radii [7]. While the Crab pulsar has a magnetic field strength at the 
light cylinder of 9.3 xlO 5 G [6], this value ranges from only 3 to 30G for these sources, 
suggesting that the bursts originate from a different emission mechanism. They therefore 
concluded that these sources represent a previously unknown population of bursting neutron 
stars, which they call rotating radio transients (RRATs). These interesting new observa- 
tions make us to investigate the rotation (e.g., it is still controversial how much angular 
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momentum the iron cores have before the onset of the gravitational-collapse) [1,6] as well 
as Pauli-blocking effects [8] in compact stars considering the phase diagram [4]. 
In present approach the Uehling-Uhlenbeck collision term [9-10] which could describe the 
collision of a gas of dilute hard-sphere Fermi- or Bose-particles by tuning a parameter 7 : 
a Pauli-blocking factor (or 7/ with / being a normalized (continuous) distribution function 
giving the number of particles per cell) is adopted together with a free-orientation 9 (which 
is related to the relative direction of scattering of particles w.r.t. to the normal of the prop- 
agating plane-wave front) into the quantum discrete kinetic model [10] which can be used 
to obtain dispersion relations of plane (sound) waves propagating in different-statistic gases 
(of particles). We then study the critical behavior based on the acoustical analog [11-12] 
which has been verified before. The possible phase diagram (as the orientation is changed) 
we obtained resemble qualitatively those proposed before [4]. 



II. THEORETICAL FORMULATIONS 



The velocities of particles are restricted to, e.g., u 1; u 2 , • • • , u p , p is a finite positive 
integer. The discrete number densities of particles are denoted by ATj(x, t) associated with 
the velocities Uj at point x and time t. If only nonlinear binary collisions and the evolution 
of iVj are considered, we have 

^ + n i -VN i ^F i = W{A%N k N l -Af j N i N j ), i e A = {1, ■ • • ,p}, (1) 

j,k,l 



where and (k,l) (i 7^ j or k 7^ I) are admissible sets of collisions 10]. Here, the 

summation is taken over all j, k, I e A, where A 1 ^ are nonnegative constants satisfying Q] 
Am = Al x = Ajl, + u, — u& — uj) = 0, and A% = A®. The conditions defined 

for the discrete velocities above require that there are elastic, binary collisions, such that 
momentum and energy are preserved, i.e., u, + Uj = + u;, |iij| 2 + |u.,-| 2 = |u^| 2 + |u^| 2 , are 
possible for 1 < k,l < p. We note that, the summation of iV, Ni) : the total discrete 
number density here is related to the macroscopic density : p (= m p J2i Ni), where m p is the 
mass of the particle 110]. 

n 

Together with the introducing of the Uehling-Uhlenbeck collision term |10J] : F{ = 
Ej, k ,l A% [NM (1 + 7 JVi)(l + 7^) -NiN.il + 7 jV fc )(l + 7 JVj)], into equation (1), for 7 < 
(normally, 7 = —1), we can then obtain a quantum discrete kinetic equation for a gas of 
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Fermi-particles; while for 7 > (normally, 7 = 1) we obtain one for a gas of Bose-particles, 
and for 7 = we recover the equation (1). 

Considering binary collisions only, from equation above, the model of quantum discrete ki- 
netic equation for Fermi or Bose gases proposed before is then a system of 2n(= p) semilinear 
partial differential equations of the hyperbolic type : 
d d cS 2n 

di Nt + Vl ' te Ni = ~ E^'^W 1 + 7^+0(1 + 7^+n+l)" 

2cSN i N i+n (l + 1 N i+1 )(l + 1 N i+n+1 ), i = l,---,2n, (2) 

where iVj = Ni +2n are unknown functions, and Vj =c(cos[# + (i — l)n / n], sin[9 + (i — 1)tt / n}) ; 

c is a reference velocity modulus and the same order of magnitude as that used in Ref. 10 

(c, the sound speed in the absence of scatters), 9 is the orientation starting from the positive 

x— axis to the u\ direction, S is an effective collision cross-section for the collision system. 

Since passage of the plane (sound wave) will cause a small departure from an equilibrium 

state and result in energy loss owing to internal friction and heat conduction, we linearize 

above equations around a uniform equilibrium state (particles' number density : No) by 

setting Ni(t,x) —Nq (1 + Pi(t,x)), where Pi is a small perturbation. After some similar 

manipulations (please refer to Chu in [12-13]), with B = ^yN < 0, which gives or defines 

the (proportional) contribution from the Fermi gases (if 7 < 0, e.g., 7 = —1) or the Bose 

gases (B > 0, if 7 > 0, e.g., 7 = 1), we then have 

r <9 2 2 2r „ (m-lW.d 2 . , n J ln 4cSN (l + B) " d „ 

+ c i ox > [e + + icSNa(l + B) -]D m = °± > g -ft, (3) 

where D m = (P m + P m+n )/2, m = 1, • • • , n, since Di = D m for 1 = m (mod 2n). 
We are ready to look for the solutions in the form of plane wave D m = a m exp i(kx — out), 
(m = 1, • • • , n), with uj—uj{k). This is related to the dispersion relations of ID (forced) plane 
wave propagation in Fermi or Bose gases. So we have 

(1 + ih(l + B)- 2\ 2 cos 2 [6 + (m ~ 1)?r ]K - ^ + B) E a k = 0, m = 1, • • • , n, (4) 

n n k=1 

where 

A = kc/(V2uj), h = AcSNq/uj oc l/K n , 

where h is the rarefaction parameter of the gas; K n is the Knudsen number which is de- 
fined as the ratio of the mean free path of gases to the wave length of the plane (sound) wave. 
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III. NUMERICAL RESULTS AND DISCUSSIONS 



We firstly introduce the concept of acoustical analog [11-12] in brief. In a mesoscopic 
system, where the sample size is smaller than the mean free path for an elastic scattering, it 
is satisfactory for a one-electron model to solve the time-independent Schrodinger equation 



: -(h 2 /2m)V 2 tp + V'(r)ip = Etp or (after dividing by -Ti 2 /2m) V 2 ^ + [q 2 - V{f)]ip = 0, 
where q is an (energy) eigenvalue parameter, which for the quantum-mechanic system is 



or (after applying a Fourier transform in time and contriving a system where c (the wave 
speed) varies with position r) V 2 tp + [q 2 — V(r)]t/j = 0, here, the eigenvalue parameter q is 
cu/co, where oj is a natural frequency and Co is a reference wave speed. Comparing the time 
dependencies one gets the quantum and classical relation E = hu. The localized state could 
thus be determined via E or the rarefaction parameter (h) which is related to the ratio of 
the collision frequency and the wave frequency [11-12]. 

The complex spectra (A = A r + i A*; the real part A r = k r c/(\/2uS): sound dispersion, 
a relative measure of the sound or phase speed; the imaginary part A« = kic/(y/2uj) : 
sound attenuation or absorption) could be obtained from the complex polynomial equation 
above. Here, the Pauli-blocking parameter (B) could be related to the occupation number 
of different-statistic particles of gases [10]. To examine the critical region possibly tuned by 
the Pauli-blocking measure B = jNq and the free orientation 0, as evidenced from previous 
Boltzmann results [12] : A« = for cases of 9 = it/ '4 (or B = —1), we firstly check those 
spectra near = 0, say, = 0.005 and = 7r/4 ps 0.7854, say, = 0.78535 for a 5-sweep 
(B decreases from 1 to -1), respectively. Note that, as the free-orientation is not zero, 
there will be two kinds of propagation of the disturbance wave : sound and diffusion modes 
[13-14]. The latter (anomalous) mode has been reported in Boltzmann gases (cf. Ref. 12 by 
Chu) and is related to the propagation of entropy wave which is not used in the acoustical 
analog here. The absence of (further) diffusion (or maximum absorption) for the sound 
mode at certain state (h, corresponding to the inverse of energy E; cf. Chu in Ref. 12) is 
classified as a localized state (resonance occurs) based on the acoustical analog [12]. The 
state of decreasing h might, in one case [15], correspond to that of T (absolute temperature) 
decreasing as the mean free path is increasing (density or pressure decreasing). 
We have observed the max. A« (absorption of sound mode, relevant to the localization 




2 . Meanwhile, the equation for classical (scalar) waves is V 2 ?/> — {d 2 i\)/c 



■ 2 dt 2 ) = 
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length according to the acoustical analog [12]) drop to around four orders of magnitude 
from 9 = 0.005 to 0.78535 (please see Chu (2001) in [12])! This is a clear demonstration 
of the effect of free orientations. Meanwhile, once the Pauli-blocking measure (B) increases 
or decreases from zero (Boltzmann gases), the latter (Fermi gases : B < 0) shows opposite 
trend compared to that of the former (Bose gases : B > 0) considering the shift of the max. 
Aj state : 5h. 5h > is for Fermi gases (\B\ increasing), and the reverse (5h < 0) is for Bose 
gases (B increasing)! This illustrates partly the interaction effect (through the Pauli exclu- 
sion principle). These results will be crucial for further obtaining the phase diagram (as the 
density or temperature is changed) tuned by both the free orientation and the interaction. 
Here, B = — 1 or # = 0, 7r/4 might be fixed points. 

To check what happens when the temperature is decreased (or h is decreased) to near T = 
or T = T c , we collect all the data based on the acoustical analog from the dispersion relations 
(especially the absorption of sound mode) we calculated for ranges in different degrees of the 
orientation (here, 9 is up to tt/4 considering single-particle scattering and binary collisions; 
in fact, effects of 9 are symmetric w.r.t. 9 = vr/4 for < 9 < tt/2; cf. Chu (2001) in Ref. 12) 
and Pauli-blocking measure. After that, we plot the possible phase diagram for the inverse 
of the rarefaction parameter vs. the orientation (which is related to the scattering) into Fig. 
1 (for different Bs : B = -0.98,-0.9,-0.1,0,0.1,1). Here, the Knudsen number (K n ) oc 
MFP/A S with MFP and \ s being the mean free path and wave length, respectively and the 
temperature vs. MFP relations could be, ine one case, traced from Ref. 15 (following Fig. 
3 therein). This figure shows that as the temperature decreases to a rather low value, the 
orientation will decrease sharply (at least for either Bose or Fermi gases). There is no doubt 
that this result resembles qualitatively that proposed before [16] for cases of the pressure vs. 
temperature analogue. 

On the other hand, qualitatively similar results (cf. Fig. 4 in [17] : theirin higher relative 
pressure corresponding to A r here) show that (i) once the orientation 9 increases, for the same 
h (or temperature [15]), the dispersion A r (or the relative pressure [17]) increases (please 
refer to Chu (2001) in [12]); (ii) as \B\ (B: the Pauli-blocking parameter) increases, the dis- 
persion (A r ) will reach the continuum or hydrodynamical limit (larger h or high temperature 
regime) earlier. The phase speed of the plane (sound) wave in Bose gases (even for small 
but fixed h) increases more rapid than that of Fermi gases (w.r.t. to the higher temperature 
conditions : larger h) as the relevant parameter B increases. For all the rarefaction measure 
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(h), perturbed plane waves propagate faster in Bose-particle gases than Boltzmann-particle 
and Fermi-particle gases (e.g., see [18] or [13]). In fact, the real part (A r ) also resembles 
qualitatively those reported in [19] for T > T c cases. 

As for the imaginary part (A«), there is the maximum absorption (or attenuation) for certain 
h (the rarefaction parameter) which resembles that reported in [20] (cf. Fig. 1 (b) therein). 
We observed a jump of the (relative) sound speed in the multiple scattering case (cf. Chu 
(2002) in [12]) which was also reported in [21] (at the phase transition between hadronic 
phase and QGP). 

To know the detailed effects of interactions (tuned by the Pauli-blo eking parameter : B here) 
and the orientation, which could be linked to the effective number of thermodynamic degrees 
of freedom (as stated in [22] : u, for an ideal gas of massless, non-interacting constituents, 
v counts the number of bosonic degrees of freedom plus the number of fermionic degrees 
of freedom weighted by 7/8), we plot 9 (of which the localization or resonance occurs for 
specific B) vs. /i (the chemical potential, in arbitrary units) in Figs. 2 and 3 by referring 
to two possibly localized states (9 = or 7r/4; cf. Chu (2001) in [12]). The trends here 
resemble that of Fig. 1 in [16]. The density rises from the onset of nuclear matter through 
the transition to quark matter as illustrated in Fig. 2 (cf. the left branch of Fig. 1 in [16] by 
Alford). The compact star is possibly in this region of the phase diagram. Then there might 
be different behaviors separated by crossover regions as shown in Figs. 2 (9 = dominated, 
T (the temperature) oc h (the rarefaction measure)) and 3 (9 = it/ '4 dominated, T (the 
temperature) oc l/h). We remind the readers that for the case of 9 = ir/4 dominated (Fig. 
3), as the lower temperature is associated with the higher density, fermions (B < 0) link to 
the lower temperature regime (under the same orientation). 

If there are differences between ours and those reported before (presumed that the mat- 
ter under study is in (approximate) local thermal equilibrium. At RHIC, such evidence is 
believed to be provided by the agreement of the elliptic flow [23] measured in noncentral col- 
lisions with hydrodynamic model predictions. Such predictions are based on the assumption 
that the matter behaves like a fluid in local thermal equilibrium, with arbitrarily short mean 
free paths and correspondingly strong interactions), one possible explanation could be that 
the assumption of a completely thermal medium is a simplification [24]. The acoustic per- 
turbations we treated are close to the thermodynamic equilibrium (for Bose or Fermi gases). 
Other reasoning is related to the different types of particles (with or without fragmentation) 
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being considered. One interesting observation is that the attenuation of jet (quenching) 
observed at RHIC resembles qualitatively the attenuation of plane (sound) waves (cf. Chu 
in [12]). 

With above results, then our approach could provide an effective theory based on the op- 
posite picture of very strong interaction (via the tuning of B and 9) and very small mean 
free paths (h is large) [24]. This can also be useful to the study of problems in astrophysics 
: like compact stars. For instance, one of the most striking features of QCD is asymptotic 
freedom: the force between quarks becomes arbitrarily weak as the characteristic momentum 
scale of their interaction grows larger. This immediately suggests that at sufficiently high 
densities and low temperatures (corresponding to the case of Fig. 3 here; cf. Chu (2001) in 
[12] since = tt/4 is also possible and thus dominates the localized behavior or transition) 
matter will consist of a Fermi sea of essentially free quarks, whose behavior is dominated by 
the freest of them all: the high-momentum quarks that live at the Fermi surface. 
To conclude in brief, our illustrations here, although are based on the acoustical analog of 
our quantum discrete kinetic calculations, can indeed show the Fermi and Bose liquid (say, 
Cooper pairs) and their critical behavior for the transition (at least valid to the regime 
T > T c considering the QGP) once the orientation is tuned as well as the temperature is 
decreased significantly. We shall investigate more complicated problems in the future [25-27] 
(e.g., the saturated orientation shown in Fig. 3 which is almost the same for all different- 
statistic gases of particles might be relevant to the Stefan-Boltzmann limit; cf. Fig. 5 in 
[26] by Rischke). 
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Fig. 1 Possible phase diagram for different-statistic gases w.r.t. the free 
orientation (9) and Knudsen number (K n oc the mean free path/wave length). 
B > : bosonic particles; B < : fermionic particles [10]. The orientation is 
related to the scattering. K n might be transformed to the dimensionless or relative 
temperature (cf. Fig. 3 in [15] ). 



12 



10 1, 



10"' 

o 

■£10 

O 

10" 3 



10" 



Pauli-Blocking parameter 



<> ^ 


* V 


o 


B = 






X 


B = 0.1 





% v 


V 


B = 1 




o 


□ 


B = - 0.1 


A 


°5 v 





B = -0.9 




o 


A 


B = - 0.98 



OX V 



e = : dominated 



10" 



1/h oc jj, 



10* 



Fig. 2 Possible phase diagram for different-statistic gases w.r.t. the 
orientation (9) and temperature T oc h (the rarefaction measure). 
The unit of \i (related to the chemical potential) is arbitrary. B > : bosonic particles; 
B < : fermionic particles [10]. 
The trend here resembles the left branch of Fig. 1 in [16] by Alford. 
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Fig. 3 Possible phase diagram for different-statistic gases w.r.t. the 
orientation and temperature T oc 1/h or K n (the Knudsen number). 
The unit of fi is arbitrary (cf. Fig. 3 in [15]). The trend here resembles the CFL (i 

branch of Fig. 1 in [16] by Alford. 
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